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Abstract. In the perturbative QCD with N c — > oo the amplitude for the collision of two 
heavy nuclei is expressed via dipole densities in the nuclei. Coupled equations for these densi- 
ties are derived in the configuration space. The equations are conformal invariant in absence 
of external sources. Passing to conformal basis and its possible truncation are discussed. 

1 Introduction 

In the perturbative QCD with N c — ► oo , high-energy scattering of a pointlike projectile on a 
large nucleus is described by a sum of fan diagrams constructed from BFKL pomerons and 
verteces for their splitting in two. Summation of these diagrams leads to the well-known BK 
equation [1, 2, 3], now well studied both theoretically and numerically. Generalization of 
these results to nucleus- nucleus scattering requires a symmetric treatment of projectile and 
target and obviously involves not only the vertex for the splitting of a pomeron in two but 
also fusing of two pomerons in one. Such a program was realized in our papers [4]. There 
we limited ourselves to the case when momenta transferred to both nuclei are neglegible as 
compared to gluon momenta inside pomerons, which physically corresponds to the limit of 
very heavy nuclei and a finite pomeron slope. As a result all pomerons were propagating 
in the forward direction. This simplified the final equations considerably but their basic 
conformal invariance property remained hidden. Also written in the momentum space the 
equations are difficult to compare with the results following from the dipole picture and so- 
called JIMWLK approach, in which the problem of symmetric treatment of projectile and 
target (and also of inclusion of pomeron loops) is lately being studied very actively (see e.g. 
[5] and references therein). For these reasons in this paper we rederive equations describing 
nucleus-nucleus scattering for a general case, when the pomerons are allowed to change their 
momenta in their interaction with the nuclei and between themselves. Our final equations 
are in the transverse coordinate space, so that comparison with the dipole approach will be 
facilitated. 

It is important to stress the approximations used in the derivation. We rely on the 
perturbative QCD in the limit of large number of colours, iV c — > oo, and assume both nuclei 
to be large, with their atomic numbers A,B >> 1. Thich allows to take into account only 
tree diagrams constructed of pomerons and their interaction verteces (see Fig. 1 for the 
simplest examples). Both pomeron loops and contributions which cannot be described in 
terms of pomerons (e.g. gluonic interaction between two pomerons) are neglected, since they 
are damped by higher powers of 1/N C and/or 1/A, 1/B. Obviously this approximation cannot 
work at superlarge energies when the pomeron can propagate to transverse distances larger 
than the nuclei dimensions, which occurs at rapidities Y such that a'Y > R\g, where a' 



Figure 1: Simplest tree diagrams for nucleus-nucleus scattering 



is the pomeron slope. All experimental data give a < 0.01 fin 2 , so that these rapidities 
seem to be well beyond our present and predictable possibilities. The approximation of large 
nuclei also allows to neglect correlations between colour distributions of different nucleons 
and excludes diagrams in which a pomeron interacts simultaneously with two nucleons of the 
target or projectile. 

Our final result is a pair of equations for the dipole densities in the two nuclei, which 
possess full conformal invariance in absence of external sources. The equations include terms 
with the interaction between the two densities. If one neglects this interaction the equations 
decouple into a pair of BK equations for the projectile and target. The equations are in fact 
very complicated and not well suited for numerical studies, which are difficult already for 
the forward propagation case considered in [4]. Some simplification seems to be possible 
by passing to the conformal basis and restricting to lowest conformal states, which is also 
discussed in the paper. 

2 Pomeron diagrams and effective field theory 

At fixed overall impact parameter b the AB amplitude A(Y, b) can be presented as an expo- 
nential of its connected part: 



The dimensionless T is an integral over two impact parameters 6 a and bs of the collision 
point relative to the centers of the nuclei A and B: 



As mentioned in the Introduction, in the perturbative QCD with N c — > oo the amplitude 
—T(Y, bA, bs) is given by a sum of all connected tree diagrams constructed of BFKL pomerons 
and the triple pomeron vertex. More concretely, in these diagrams a line ("propagator") 
describing propagation of a pair of gluons from rapidity y and points r\ and r 2 to rapidity y' 
and points r[ and r' 2 corresponds to the BFKL Green function G y - y '(ri, r2\r[, r' 2 ) [6]: 
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Figure 2: Splitting of a pomeron into two (a) and fusion of two pomerons into one (b). 

where fi = {n, v, ro} = {h, ro}, summation in (3) includes summation over n and integrations 
over v and transverse vector ro with the weight {y 2 + n 2 /4)/ir 4 ; functions form the 
conformal basis. In the complex notation r = x + iy, r* = x — iy 



rig y 



12 



10' 20 



(4) 



where r\ 2 = r\ — r 2 etc; h = (1 — n)/2 + iv and h = 1 — h* are conformal weights. Function 
uj^ is the BFKL eigenvalue 



n\ + 1 

Y>(1) - Reip[ J — ^ h ii> 



(5) 



where standardly a = a s N c /ir. Finally 



A» = Aft 



[(n + 1) 2 + 4i/ 2 ][(n- l) 2 + 4z/ 2 ]" 



(6) 



The interaction between pomerons is realized via the triple pomeron vertex. It is non- 
local and not symmetric in the incoming and outgoing pomerons. In the limit N c — > oo its 
form for the splitting of a pomeron into two was established in [7]. The three BFKL Green 
functions are connected by it as follows (see Fig. 2a) 



2a 2 s N c f d 2 nd 2 r 2 d 2 r 3 

lyi 2 yi 2 /I-, 2 

'12' 32' 13 



KNc r 
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Here L\ 3 is (up to a numerical factor) the quadratic Kasimir operator for the conformal 
transformations of r\ and r 3 : 

Lis = r\ zV \vl (8) 

where in the configuration space p 2 = — V 2 . In (8) l\ 3 is acting on the left. Note that the 
triple pomeron vertex is symmetric in the gluons inside the outgoing pomerons (i.e under 
r\ <-> r 2 and r 3 <-> r 2 ). So the outgoing pomerons have to be symmetric in their respective 
gluons. The form of the vertex for the fusion of two pomerons into one is actually not known. 
However, the symmetry between target and projectile prompts us to assume that for the 
inverse process 2 + 3^1 the BFKL functions are to be joined as (Fig. 2b) 

1o? a N c f d 2 r x d 2 r 2 d 2 r 3 , , „ „ , , 

— — / 2 2 2 G y' 2 -y( r i,r 2 \ri,r2)G y , y(r 3 ,r 2 \r 3 ,r 2 )L 13 G y _ y ,{r 1 ,r 3 \r l ,r 3 ). (9) 

n J r 12 r 32 r 13 

Finally we have to describe the interaction of the pomerons with the two nuclei. The 
BFKL Green functions corresponding to the external legs of the diagrams are to be integrated 
with the colour density of each nucleus. We take the target nucleus (B) at rest, that is, at 
rapidity zero. Then each outgoing external BFKL Green function is to be transformed into 

J d^d^Gyin^yjpBir'^) = J dtf#r> 1 (Pi J 2 G v ^{r u r2\i J 1 S 2 )TBh/AA), (10) 

where p B is the colour density of the target. If we neglect correlations between the colours 
of the nucleons in the nucleus B 

r B (y,ri,r 2 ) = g 2 BT B (b B )5 2 (b B - H±Z^ p N (r 12 )5(y), (11) 

where T B is the profile function of the nucleus B and pat is the colour density of the nucleon. 
Similarly each ingoing BFKL external Green function is to be transformed into 

J dy'd 2 r[d 2 r' 2 T A {y'A,r2)G y ,- y {r[y 2 \n,r 2 ), (12) 

where 

T A (y,n,r 2 ) = g 2 AT A (b A )5 2 (b A - H±I^j p N (r 12 )5(y - Y) (13) 

and Y is the overall rapidity. If the densities t a and t b are symmetric in the gluons (which 
is true for (11) and (13) and will be assumed in the following) then the initial pomerons are 
also symmetric in the gluons inside them. Together with the mentioned properties of the 
triple pomeron verteces it means that only symmetric pomeron states are propagating in the 
two nuclei. 

To find the amplitude, one has to sum over all connected diagrams with M ingoing and N 
outgoing lines, corresponding to M interactions with the projectile and N interactions with 
the target, divided by MINI. 

It is trivial to see that this sum exactly corresponds to the sum of tree diagrams generated 
by an effective quantum theory of two pomeronic fields (dipole densities, up to a factor) 
&(y,ri,r 2 ) and &(y, r±, r 2 ), symmetric in r±,r 2 , with the action 

S = So + Si + Se 

consisting of three terms, which correspond to free pomerons, their mutual interaction and 
their interaction with external sourses (nuclei) respectively. 
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To give the correct propagators So has to be chosen as 

So = J W^^dV^M^t^^!,^)^^^,^!^,^)^^,/!,^) EE (^IG- 1 ^), 

(14) 

where (|) means the integration over y and both gluon coordinates. Note that the sign of So 
corresponds to the following substitution of the conventionally defined field variables <E> and 

(15) 

which allows to make all terms of the action real. 

According to (7) and (9) the interaction term Si is local in rapidity 



gj= 2o#Ve f dy f cPnd^cPrz | ,^ r . ^ ^ ra)$f ( ^ ^ + ^ c | (lg) 
ir J J r 12 r 32 r 13 >- \ 



The overall sign combines the initial factor z and i 3 from the substitution (15). 
Finally the interaction with the nuclei is local both in rapidity and coordinates: 

Se = ~ J dyd 2 r 1 d 2 r 2 {$(y,r 1 ,r 2 )T A (y,ri,r 2 ) + &(y, n, r 2 )r B (y, n, r 2 )}. (17) 

The minus sign comes from the initial % and the substitution (15). 

The amplitude T(Y, &b) is then expressed through a functional integral 

Z = [ D$D&e s . (18) 



In the tree approximation, corresponding to diagrams of the type shown in Fig. 1, keeping 
only connected diagrams we find 

T(Y, b A , b b ) = -ln^- = $t} (19) 

where Zq{= 1) is the value of Z for Se = 0. and the action 5 is to be calculated for $ and 
satifying the classical equation of motion. 

3 Equations for the dipole densities 

The classical equations of motion follow from the variation of the action with respect to <E> 
and 

5S _ 5S _ () 
<S$(y,ri,r 2 ) ^(y, n, r 2 ) 

We find a pair of equations 

G~ 1 ^(y,ri,r 2 ) + 2asiVc / 2 d 2 r3 2 ri,r 3 )$(y,r 2 ,r 3 )Li 2 

+2&(y,r 3 ,r 2 )L 13 $(y,r 1 ,r 3 )} = T B (y, ri ,r 2 ) (21) 

$ t d/,ri,r 2 )G- 1 + ^^ / 2 d J" 3 2 n, r 3 )$t ( y , r 2 , r 3 )L i2 

71 ^ r 12 r 32 r 13 

+2$(y,r 3 ,r 2 )Li3$ t (y,ri,r 3 )} = T A (y,n,r 2 ), (22) 



and 
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where the operators L\ 2 are assumed to act on the left. These equations have also to be 
supplemented by conditions 

*(y,ri,r 2 ) = if y<0, &(y, ri ,r 2 ) = if y > Y. (23) 

To write the equations in their final form we note that the BFKL Green function G satisfies 
the equations 

^(| + ^ G =(| + Ft >^ G = 1 ' ^ 
where H is the BFKL Hamiltonian: 

H = £ ( Inp 2 + Inp 2 + \ln r\ 2 ■ pj + -1 In r\ 2 ■ p\ - 4^(1)) . (25) 

Fom this we conclude that 

G-'-A3(| + ff)-(| + H*)A3. (26) 

Multiplying Eqs. (21)and (22) by Pi 2 p 2 2 fom the left and from the right respectively we 
finally find 

fd t T \^< x 2a 2 N c fd 2 r 3 r 2 2 ^. . 

\dy J 7T J r 2 2 rf 3 

+ ^l^c L -i I ^3^2 $ t (2/)r , 3)r2)Ll3$(y)ri)r3) = TB (y, ri , r2 ) (27) 

(_|_ + tfW (y , ri , r2 ) + ^£ /" ^2$t (y)ri)r3 ) $ t (y)r2)r3) 
V vr J r$ 2 r( 3 



and 



/ -^^2$(y,r3,r 2 )L 1 3$ t (y,r 1 ,r 3 ) = r A (y, n, r 2 ). (28) 

T ^ ^32*13 

The 5-like dependence of the external sources on y together with conditions (23) imply that 
one can drop the sources in the equations and substitute them by the boundary conditions 
at y = and y = Y: 

$(y,ri,r 2 )y=o = pB(n,r 2 ), $ f (y, n, r 2 ) y=Y = pa(ti, r 2 ), (29) 

where pa,b are given by (11) and (13) without the <5-functions in rapidity. 

From the form of equations one immediately concludes that they are conformal invariant 
provided the external sources possess this invariance. 

These equations can be also written in the form which allows an easy comparison with 
the BK equation for non-forward fan diagrams. To this aim one rescales the fields putting 

$(n , r2) = ^l, , Hnin) = q^2l ( 30 ) 
4ira s 4ira s 

and uses a representation for the Hamiltonian H [8] 

#/(ri,r 2 ) = f- / ^^(/(r 1 ,r 2 )-/(r 1 ,r 3 )-/(r 2 ,r 3 )). (31) 

/7r J r 23 r 13 
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Then our equations take the form 
dN(r\,r2) a f d 2 r 3 r\ 2 



0lJ : > ;T / (N(n,r 2 ) - N( n , r 3 ) - N(r 2 , r 3 ) + iV(n, r 3 )iV(r 2 , r 3 )) 

/" 1Z^2 (r 3 ,r 2 )L 13 N( n ,r 3 ) = (32) 

7r j r^9 r ia 



and 

dN^ (r\,r 2 ) a f d 2 r 3 rf 2 



dy 2itJ r| 3 rf 3 



(iVt(n,r 2 ) -iVt(n,r 3 )- Art(r 2 ,r 3 ) + Art(n,r 3 )Art(r 2 ,r 3 )) 



+^ L -i /" !^L2 iV ( r 3, r2 )L 1 3iV t (r 1 ,r3),= (33) 

n J r 32 r 13 



with the boundary conditions which follow from (29) after rescaling (30). 

If one neglects the second term on the r.h.s in both equations and thus decouples N and 
the equations turn into a pair of independent BK equations for dipole densities in the 
nuclei A and B evolving in opposite directions in rapidity and corresponding to the sum of 
two sets of fan diagrams starting from the projectile or target. However the second terms 
introduce interaction between these sets and correspond to diagrams which contain both 
splitting and fusion of pomerons. The structure of this interaction is rather complicated in 
both configuration and momentum spaces due to non-locality of the inverse operator L^ 1 . 
One expects it to be simplified in the conformal basis, which will be the subject of the next 
section. 

Meanwhile, using the equations of motion one can simplify the expression for the action 
S calculated on their solution. Indeed multiplying Eqs. (21) and (22) by &(y, ri,r 2 ) and 
$(y,ri,r 2 ), integrating over y, ri,r 2 and summing the results one obtains a relation 

2So + 35/ + S E = 0. (34) 

This can be used to exclude one of the parts of the action when calculating the amplitude T. 
Recalling that the fields are discontinuous at the boundaries we obtain from (34) 

T(Y, b A , b B ) = ^(S E - S ) = l(Si- S E ). (35) 



4 Equations in the conformal basis 

One may hope that the equations for the gluon densities may be somewhat simplified in the 
conformal basis formed by functions £^ M (ri,r 2 ). To this end we present 

Hy,r 1 ,r 2 ) = J2E fl (r 1 ,r 2 )^^y). (36) 
The orthonormalization properties of £ , M (ri,r 2 ) [6] allow to invert this relation and find 

%(y) = I ^^^(ri,r 2 )*(j/,n,r 2 ). (37) 

J r 12 

Since [i = {n,is, ro}, transition to the conformal basis by itself does not change the number 
of variables (three). However it drastically simplifies the operators L in the the mixing term 
of our equations. 
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Indeed the mixing term of Eq. (27) can be written as 

T mix (n,r 2 ) = ^±L- X l [ £ d>t i ( y )d, M2( y)A M2 1 ^ i (r3,r 2 )^ 2 (r 1 ,r3), (38) 

17 J r 32 r 13 M1 , M2 

where we have used that 

L 13 E^( ri ,r 3 ) = X^E^ ri ,r 3 ). (39) 

Expanding the integral over r 3 considered as a function of r\ and r 2 in the conformal basis 
we get 

T mix (n,r 2 ) = ^T™^( n ,r 2 ), (40) 

where according to (39) 

4a 2 N c f d 2 rid 2 r 2 



L i2 I -S^ E ®lM*MKt E lMir2)E^{ ri ,n). (41) 

•/ '39!' 13 ... .... 



'32' 13 miAt2 

We integrate by parts transforming action of L^ 1 on E^(ri,r 2 )/rf 2 and use 

rf 2 L^r^ = L- X l (42) 
to apply directly on £^(ri,r 2 ) which gives a factor A M . So in the end we get 
Aa 2 N c f d 2 r 1 d 2 r 2 d 2 r 3 

"2 >. a 2 j. a 2 

12 r 32 r 13 ^ 1)/12 

4«2Ar c 



T T = ^f d2 lffi n ^ E >^) E <(y)^ 2 (y)A; 2 1 ^ 1 (r 3 ,r 2 )^ 2 (r 1 ,r 3 ) 

•/ ' 19' 32' 13 „, 



a m E ^.w^lfo)*^). (43) 



7T 



where V^ )M1)/12 is the triple pomeron vertex in the conformal basis 

/d^Txffiv d 2 T 3 
g-2-2 — Epin, r 2 )E^ (r 3 , r 2 )^ 2 (n, r 3 ) (44) 
r 12 r 32 r 13 

and /x corresponds to the complex conjugate basis function: if fx = {h, ro} then fl = {1 — h, ro} 
(and always h = 1 — h*). 

A similar transformation of the first integral term in Eq. (27) is straightforward and leads 
to the result which is different from (43) by the absence of conjugate $'s and factors A. So 
we find the first equation in the conformal basis as (suppressing the common argument y and 
dropping the source term) 

=-A"^ E *4W 2 ^+2^W 2 <) (45) 

°y 71 hum A ^ 

The second equation can be obtained by reversing the direction of propagation in rapidity 
and passing to conjugate fields: 

= -^ + ^£ £ 4,(^^ + 2^^^). (46) 
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The triple pomeron vertex V^ ;jU1)/12 was studied in [9]. It depends on three conformal 
weights h, hi and /i 2 and three center-of-mass vectors {ro, roi, ro 2 } = {po 5 Pi 5 P2j-- The depen- 
dence on the latters is determined by the conformal invariance, so that (in complex notation) 

v m2 = n h , hlM n ,>;, A ■ (47) 

i<j 

where i = 0, 1,2, A01 = ho + hi — /i 2 , A01 = ^0 + hi — J12 etc. The part of the vertex 
depending on conformal weights QhMM was found in [9] for arbitrary conformal weights 
in terms of the Meijer function G44. The complicated form of the vertex together with the 
use of complex variables make equations (45), (46) for the gluon densities in the general case 
not very suitable for practical calculations, in spite of the simplification for the action of 
operators L. However they may serve as a starting point for further simplifications realized 
by truncating the equations by certain low values of conformal weights. In the next section 
we consider a most drastic example of such a truncation. 



5 Lowest conformal weights 

As well-known from the study of the linear BFKL equation in the high-energy limit, the 
leading contribution comes from the minimal conformal weight in the expansion (36) , namely 
h = h = 1/2, which corresponds to n = and v = 0. So the simplest case which may be of 
interest for our problem is to put n = in all places and v = whenever this is allowed by 
the equations, that is in <E>^, <E>t and £lh h^- Then one finds for the triple pomeron coupling 
[9] 

^1/2,1/2,1/2 = Vo = 2^ 7 4 F 3 (l/2) 6 F 5 (l/2) = 7766.679. (48) 

The unknown fields and ^ become functions of rapidity y and center-of mass vector po 
(actually of p"^ due to rotational invariance). The two equations (45) and (46) simplify to 

d*(f>o) = 
dy 

a2 s N c^ f d 2 pid 2 p 2 poi 



A*(P0) - ^0 / lM^^(i n _^_)^(i n _5W_)< &(pi )($( p2 ) + 2 d>t (p2) ) (49) 

OTT' J P01P02 V P02P12 7 V P01P12 7 V ' 



and 



P01P02 v P02/?12 / v P01P12- 



-A*t(p ) + ^o / ^^^(ln-^_)^(ln-^-)*t (pi )(*t ( ^ ) + 2 * (p2) ) > (50) 

8lT< J P01P02 V P02P12 7 V P01P12 7 v ' 



dy 

a 2 s N Cr ^ f d 2 pid 2 p 2 s „ 1 ln poi 

P01P02 ^ PmPvi' v PoiPw 
where A = uj n= o tU= o is the BFKL intercept. We have taken into account that due to the 
presence of the 8- functions we have P12 = 1 (in the chosen scale, determined by the sources). 

From the assumed independence of the fields of v it follows that the boundary conditions 
for these equations have to belong to the class of functions of the form 



/(ri,r 2 )=/^y'(ln^W 3) , (51) 

where g(r) is an arbitrary function. Obviously this restricts the sources to be of a very special 
sort, with the dependence on two vectors r\ and r 2 and thus on three variables r 2 , r\ and rir2 
determined by a function of a single variable r 2 . So any practical use of the thus simplified 
system of equations is questionable. At most it may serve to study the qualitative features of 
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the solution in the limit of high energies, when one may hope that the influence of the choice 
of the boundary conditions becomes neglegible (as it happens with the BK equation). Still 
even Eqs. (49) and (50) do not look easily solvable. We reserve their study for a separate 
publication. 

6 Conclusions 

We have derived a pair of equations for dipole densities in two heavy nuclei, which describe 
their scattering in the perturbative QCD with a large number of colours. The equations 
contain mixing terms which are both non-linear and non-local. In absence of mixing the 
equations decouple into a pair of BK equations for the projectile and target. 

In contrast to the hA case the equations are to be solved with given boundary conditions 
at rapidities of the projectile and target, which complicates their solution enormously. The 
equations themselves are conformal invariant. This invariance is naturally broken by the 
sources. However use of the conformal basis may open ways for various simplifications of the 
equations, which may facilitate their solution, if only on the qualitative level. 
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